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ABSTRACT 

The spectral variation of the cosmic microwave background (CMB) as ob- 
served by WMAP was tested using foreground reduced WMAP5 data, by producing 
subtraction maps at the 1° angular resolution between the two cosmological bands 
of V and W, for masked sky areas that avoid the Galactic disk. The resulting 
V — W map revealed a non-acoustic signal over and above the WMAP5 pixel noise, 
with two main properties. Firstly, it possesses quadrupole power at the ~ 1 fiK 
level which may be attributed to foreground residuals. Second, it fluctuates also 
at all values of £ > 2, especially on the 1° scale (200 ^ i ^ 300). The behavior is 
random and symmetrical about zero temperature with a r.m.s. amplitude of ~ 
7 fiK, or 10 % of the maximum CMB anisotropy, which would require a 'cosmic 
conspiracy' among the foreground components if it is a consequence of their ex- 
istences. Both anomalies must be properly diagnosed and corrected if 'precision 
cosmology' is the claim. The second anomaly is, however, more interesting be- 
cause it opens the question on whether the CMB anisotropy genuinely represents 
primordial density seeds. 

1. Introduction 

Studies of the cosmic microwave background (CMB, Penzias and Wilson (1965)), the 
afterglow radiation of the Big Bang, are currently in a period of renaissance after the break- 
through discovery of anisotropy by the COBE mission (Smoot et al (1992)). Confirmed with 
much improved resolution and statistics by WMAP (Hinshaw et al (2009)), the phenomenon 
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provides vital information on the primordial 'seeds' of structure formation. The anisotropy 
is attributed to frequency shift of CMB light induced by these 'seed' density perturbations, 
which has the unique property that it leads to changes in the temperature of the black body 
spectrum and not the shape of it. The CMB has maximum anisotropy power at the 1° scale, 
or harmonic number ^ « 220, with lower amplitude secondary and tertiary peaks at higher 



The ACDM cosmological model (Spergel et al(2007)) explains the entire power spectrum 
remarkably using six parameters, by attributing the peaks to acoustic oscillations of baryon 
and dark matter fluids, as long wavelength modes of density contrast enter the horizon and 
undergo causal physical evolution. CMB light emitted from within an overdense region of 
the oscillation are redshifted by a constant fractional amount, resulting in a cold spot, which 
is a lowering by 5T of the black body temperature T, and is frequency independent, i.e. 
5T/T — bvjv — constant. The opposite effect of blueshift apphes to underdense regions, 
leading to hot spots. Therefore, if the anisotropy is genuinely due to acoustic oscillations, 
the inferred change in T at a given region should be the same for all the 'clean' frequency 
passbands of the WMAP mission. Since a corresponding variation of the CMB flux B(v^ T) at 
any given frequency u is 6B = {dB/dT)6T if the cause is solely 6T with no accompanying 
distortion of the functional form oiB itself, the expected 6B at constant 6T is then the 'dipole 
spectrum' dB/dT which is well measured by COBE-FIRAS (Mather et al(1994)). Moreover, 
the WMAP data are cahbrated w.r.t. this dipole response. 

A noteworthy point about the acoustic peaks is that one needs to employ the technique 
of cross correlation to reduce the noise contamination at high £, especially the harmonics 
of the second and higher acoustic peaks. Specifically one computes the all-sky cross power 
spectrum 



where the indices i and j denote independent data streams with uncorrelated noise that arise 
from a pair of maps at different frequency bands (or same band but taken at different times), 
and = ST^^ is the apparent CMB temperature anisotropy for the spherical harmonics 
{i,m) as recorded by observation i. Since the use of multiple passbands is crucial to the 
accurate profiling of the acoustic oscillations, it is important that we do compare them with 
care, down to the level of measurement uncertainties. Only a priori statistically consistent 
maps should be cross correlated, in the sense that any real discrepancies between such maps 
may carry vital information about new physical processes that their cross power spectrum 
does not reveal. In one previous attempt to address this point (see Figure 9 of Bennett et 
al(2003a)) WMAPl data downgraded to an angular resolution commensurate with COBE were 
used to produce a difference (subtraction) map between the two missions. When displayed 




(1) 



m 



- 3 - 



side by side with the map of the expected noise for each resuhing pixel, the two maps did 
appear consistent. Nevertheless, this powerful method of probing the CMB anisotropy does, 
in the context of the specific datasets used by Bennett et al(2003a), suffer from one setback: 
it is limited by the sensitivity and resolution of COBE. 

In another test of a similar kind, we observe that each amplitude a^^ can further be 
factorized as = aembl, where the array b} accounts for the smoothing effects of both the 
beam and the finite sky map pixel size, and a£m = ST^m is the true amplitude of the CMB 
anisotropy. The results (see Figure 13 of Hinshaw et al (2007)) indicate agreement of the 
variance Cf, hence ST^, within the margin of a few percent for i < 400 among the many 
cross power spectra formed by the various possible combinations of pairs of all-sky maps. 
This offers more ground for optimism, but to be definitive the remaining discrepancy needs 
to be demonstrably attributed to noise, instrumental systematics, or foreground emission. 

The purpose of our investigation is to perform further, more revealing comparisons than 
the two past ones described above, initially by focussing upon the angular scale of the first 
acoustic peak, which is ~ 1°. Our analysis will be done in both real (angular) and harmonic 
domains, because while most of the effort have hitherto been pursued in the latter, the former 
is the domain in which the raw data were acquired and organized. 

2. The all sky difference map between the WMAP5 V and W bands 

We adopted the Healpis0 pixelization scheme to ensure that all pixels across the sky 
have the same area (or solid angle). Firstly the W band data is smoothed to the V band 
resolution. Then the whole sky map is downgraded to ~ 1° diameter (corresponding to 
nside of 64 in the parametrization of the WMAP database), which is not only commensurate 
with the scale of global maximum 6T power, but also large enough to prevent data over- 
sampling due to the use of too high a resolution, as the size is comfortably bigger than the 
beam width of the WMAP V band (61 GHz) larger than that of the W band. 

The resulting 6T values for the two cosmological passbands of V and W, span f^i 35,000 
clean (i.e. ext-maskecEl and foreground subtracteclfl) pixels, from which aV — W difference 
map at this ~ 1° resolution was made. After removing the monopole and dipole residuals 
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(the latter aligned with the original COBE dipole), this map is displayed in Figure [T] along 
with the corresponding pixel noise map for reference; the latter represents the expected 
appearance of the V — W map if the CMB anisotropy is genuinely acoustic in nature, so 
that the map would consist only of null pixels should the WMAP5 instruments that acquired 
them be completely noise free. When comparing the real data map of Figure [1^ with the 
simulated map of Figure [T]d, the former appears visibly noisier on the resolution scale ~ 1°; 
moreover, the Leo and Aquarius (i.e. the first and third) sky quadrants contain more cold 
pixels than the other half of the sky, indicative of the existence of a quadrupole residual. 

The extra signals revealed by the V — W subtraction map are elucidated further in 
respect of their aforementioned properties by examining the statistical distribution of the 
pixel values across the four sky quadrants. As shown in Figure |21 the distribution of the 1° 
anisotropy is considerably wider than that expected from the WMAP5 pixel noise for all the 
quadrants, by ~ 10 fiK, which is ~ 10 % of the 75 fiK power in the first acoustic peak, 
and is therefore very significant. A detailed confirmation by Gaussian curve fitting is given 
in Table [H 

The V — W quadrupole is more subtle, and is evident in the residual plots at the 
bottom of each graph in Figure [2l from which a slight skewness of the data to the right is 
apparent in quadrants 1 and 3 (the quadrants of the CMB dipole), with 2 and 4 exhibiting 
the opposite behavior. For this reason, the effect does not manifest itself as shifts in the 
Gaussian mean value n of Table [TJ Rather, the high statistical significance of both the 
quadrupole and the degree-scale signals, with the former having a magnitude of ~ 1 nK, 
are estabhshed by computing the cross power spectra of the temperature difference maps. 
Figure m This was performed at the resolution of nside= 64 using the PolSpice softwar^. 
From Figure H] also, the presence of excess non-acoustic anisotropy at all harmonics i > 2, 
including the cosmologically important 6^1° angular scale, appears robust. At the 1° scale 
(200 ^ i ^ 300), the r.m.s. is about 7 fiK, or 10 % of the maximum CMB anisotropy. 
Lastly, the V — W quadrupole may be displayed in isolation by arranging the data of the 
subtracted map as a multipole expansion 
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and evaluating at £ = 2 the amplitude 
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(note STiiO^cj)) is always a real number if the original data ST{6,(j)) are real). The ensuing 
whole sky map is in Figure [3l and the coordinates of the axes are in Table [2l 

3. Interpretation of results 

The WMAP5 V — W map reveals two principal anomalies to be explained: (a) the 
quadrupole at ^ = 2, with an amplitude of ~ IfiK, and (b) the higher harmonic signals, 
especially the ~ 8 fiK anisotropy at £ > 200 (Figure Hj). Similar findings are also made by 
others, like the noticeable hemispherical power asymmetry in the WMAPl analysis of Eriksen 
et al (2004) and confirmed in the WMAP5 data by Hoftuft et al (2009), or the large scale 
distribution investigated by Diego et al (2009). Also because both (a) and (b) are not small 
effects, claims to precision cosmology are overstatements until they are properly accounted 
for and the cosmological model accordingly adjusted. 

Concerning (a), unlike the dipole, there is no previous known CMB quadrupole of suf- 
ficient amplitude to justify its dismissal as a cross band calibration residual. In fact, our 
reported amplitude of 1 fiK is about 7 % of the 211 /xK^ WMAP5 anisotropy in the unsubtracted 
maps of the individual bands at £ = 2, which is far larger than the calibration uncertainty 
of ~ 0.5 % (Hinshaw et al (2009)) for each band. 

It will probably be more rewarding to search for remaining foreground contamination 
not yet removed by the standard data filtering and correction procedures of the WMAP5 team 
(Bennett et al(2003b). Gold et al(2009)). Thermal dust emission might have a power law 
spectrum with an index too close to that of the Ray leigh- Jeans tail in the V and W bands 
for an appreciable V - W signal, although this is an interesting scenario worthy of further 
study (Diego et al 2009). We consider here another possibility, viz. free-free emission from 
High Velocity Clouds (HVCs, Wakker et al (2009) and references therein). The clouds are 
moving at velocities sufficiently large for any Ha emission from them to be outside the 
rangqj of the WHAM survey, the database employed to estimate the free-free contribution to 
the WMAP foreground. HVC parameters for the larger and brighter clouds can reach: ^ 
0.2 cm~^ and column density 3 x 10^^ cm~^ (Wakker et al 2008). This corresponds to an 
emission measure of two units, or 6 x 10^*^ cm~^, or ^ 0.6 fiK of V-W temperature excess 
(Finkbeiner D.P. (2003)), on par with the 1 fiK of our observed quadrupole. Moreover, as 
can be seen from the all-sky map of A^hi and an estimate of the V-W excess in Figure when 



'''Example of a HVC missed by WHAM is Hill et al 2009, a cloud of unit emission measure. A notable 
exception (counter example) would be the HVC K-complex (Haffner et al 2001), with an emission measure 
of 1.1 units, that happens to fall inside the velocity window of WHAM. 
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they are compared with Figures [3] and HI the strength and distribution of HVCs do appear 
to be responsible for a non-neghgible fraction of the observed anomaly on very large scales. 
Further work in this area is clearly necessary, and will be pursued in a future, separate paper. 

We now turn to (b), the effect that occurs on the much smaller and cosmologically most 
significant angular scale of 1°. Calibration issues are again immediately excluded here, since 
the 8 anomalous amplitude is on par with the pixel noise of WMAP5 for the scale in question 
(Tabled]). Moreover, because the subtracted V — W dipole and the (unsubtracted) V — W 
quadrupole, the latter being (a), are both relatively feeble phenomena, of amphtudes ~ 0.2 
and 1 fiK respectively as compared to the 7 fiK amplitude of (b), the prospect of smaller 
scale fluctuations having been enhanced by a larger scale one can be ruled out here. CMB 
spectral distortion during the recombination era, or subsequently from the Sunyaev-Zeldovich 
(SZ) scattering, or from other foreground re-processing that were not properly compensated 
by the data cleaning procedure of WMAP5, could all be responsible for the observed anomaly. 
Although the first two interactions (Sunyaev and Chluba (2008), Birkinshaw and Gull (1983)) 
exert much smaller influences than 7 fiK (bearing in mind that the degree of SZ needs to be 
averaged over the scale of the whole sky), the foreground could potentially play a relevant 
role in a similar way as it did at very low i. Thus, in respect of free-free emission by HVCs 
alone, until a full survey at high angular resolution is performed one cannot be certain that 
the emission measure from these clouds is too weak to account for our (b) anomaly. However, 
the action of the foreground is systematic in that it does not lead to random and symmetric 
temperature excursions (about zero) between two frequencies of V and W. More precisely, 
because the sources or sinks involved have a characteristic spectrum that differs from black 
body in a specific way, any widening in Figure [2] of the data distribution w.r.t. the expected 
simulated gaussian ought to be highly asymmetric. This obviously contradicts our findings, 
i.e. we note from Figure [2] that the widening of the data histogram is highly symmetric. As 
a result, the symptoms do not point to the foreground as responsible cause. 



4. Conclusion 

We performed a new way of testing the black body nature of the CMB degree scale 
anisotropy, by comparing the all-sky distribution of temperature difference between the 
WMAP5 cosmological bands of V and W, with their expected pixel noise behavior taken fully 
into consideration by means of simulated data. In this way a non acoustic signal is found 
in the ext-masked V — W map at the ~ 1° resolution of nside = 64, with the following 
two properties. It has a quadrupole amplitude ~ 1 fiK (Figures [21 [3l and H]) which may in 
part be attributed to unsubstracted foreground emission. It also has excess anisotropy (or 
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fluctuation) on all scales i > 2, including and especially the scales of 200 < ^ < 300 where 
most of the acoustic power resides, and about which the anomaly we reported is in the form 
of a symmetric random excursion about zero temperature with a r.m.s. ~ 8 fiK (Figures 
|5] and m Table [T]) which is ~ 10 % of the maximum acoustic amplitude found at £ ~ 220. 
This type of excursion frustrates attempts to explain the effect as foreground residuals, i.e. 
it opens the question of whether the WMAP anisotropy on the 1° scale is genuinely related to 
the seeds of structure formation. 

In any case, it is clear that both anomalies have sufficiently large magnitudes to warrant 
their diagnoses through future, further investigations, if the status of precision cosmology is 
to be reinstated. 

We are grateful to the referee for very valuable suggestions towards the improvement of 
this paper. Lyman Page, Priscilla Frisch, Gary Zank, and Barry Welsh are also acknowledged 
for helpful discussions. Some of the results were obtained by means of the HEALPix package 
(Gdrski et al (2005)). 
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V- w 




error (/xK) 


a (/iK) 


error (/iK) 


Quadrant 1 


WMAP5 


-0.23 


0.15 


16.23 


0.13 


Simulation 


0.00 


0.13 


14.70 


0.12 


Difference A 


-0.23 


0.20 


6.88 


0.40 


Quadrant 2 


WMAP5 


0.24 


0.12 


14.47 


0.10 


Simulation 


-0.04 


0.12 


12.10 


0.10 


Difference A 


0.28 


0.17 


7.94 


0.24 


Quadrant 3 


WMAP5 


-0.11 


0.16 


16.22 


0.13 


Simulation 


0.03 


0.15 


14.70 


0.12 


Difference A 


-0.14 


0.22 


6.86 


0.40 


Quadrant 4 


WMAP5 


0.40 


0.13 


14.80 


0.10 


Simulation 


-0.01 


0.13 


12.26 


0.10 


Difference A 


0.41 


0.18 


8.30 


0.23 



Table 1: Parameters for the gaussian curves that fitted the WMAP5 data and the pixel noise 
histograms (the latter are the solid lines) of Figure [2J Each parameter uncertainty is set by 
the Xmin + 1 criterion, which represents the usual 68 % (or unit standard deviation) confidence 
interval for one interesting parameter, when the error bars shown in Figure [2] are employed 
for fitting both the real and pixel noise data. The difference in the width a between the 
two models, which gives the distribution width of the additional random signal, is given by 
(A(t)^ = — erf. The smaller simulated gaussian widths for quadrants 2 and 4 (relative to 1 
and 3) is due to the higher exposure times there (which contain the heavily scanned ecliptic 
poles) leading to lower pixel noise. 



V-W quadrupole location (/, b) 


hot 


(-132.1°, -14.4°), (48.0°, 14.4°) 


cold 


(-81.5°, 68.0°), (98.5°, -68.0°) 



Table 2: Orientation of the quadrupole in the WMAP5 V-W map of Figure [3l 
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Real Data 
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Fig. 1. — The ext-masked and point sources subtracted WMAP5 V—W map, viz. the difference 
map between the CMB anisotropy as measured in the V band and the W band, for the real 
data after the removal of residual monopole and dipole components (top), and simulated 
pixel noise that reflect precisely the observational condition (bottom) . Both maps are plotted 
in Galactic coordinates with the Galactic center (/, b) = (0, 0) in the middle and Galactic 
longitude / increasing to the left. To avoid the problems of beam size variation from one 
band to the next, the W band data is smoothed to the V band resolution, then the pixels 
were downgraded to the common resolution of nside= 64 using the foreground- reduced 
WMAP5 data (see section 2); this resolution under-samples the data in both bands. The color 
scale is coded within a symmetrical range: those pixels with values beyond ±40 /iK are 
displayed in the same (extreme) color; most of such pixels are around the masked regions. 
The existence of additional non- black body signal in the real data can readily be seen from 
this comparison, as the simulated map is noticeably quieter. 
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Fig. 2. — The data points show quadrant sky occurrence frequency distribution of the differ- 
ence in the degree-scale (nside= 64) anisotropy between the WMAP5 V and W bands, while 
the errors in the data are due to the WMAP5 pixel noise for the same ext-masked quadrant 
sky area, i.e. they are the statistical fluctuations in the various parts of the solid line, which 
gives the mean histogram of this noise. The orientation of each quadrant follows the same 
convention as the sky maps of Figure [H with the 1st and 3rd quadrants marking the COBE 
dipole. 
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Quadrupole 




1.1 AiK 



Fig. 3. — V — W quadrupole of the nside= 64 WMAP5 temperature difference maps, after 
ext-masking and point source subtraction. The mathematical procedure of extracting each 
multipole i is given in eqs. (3) and (4) of the text, and the software used to do these 
computations was from anaf ast of Healpix. 
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Fig. 4. — Real and simulated (noise) power spectra of the WMAP5 V — W map. These are 
V-W cross power spectra computed by cross correlating the first three years of observations 
with the last two. The errors in the real data of the first two graphs represent the pixel noise 
power of the last graph, i.e. 4c is the average of 1,000 simulated reahzations of the V-W 
WMAP5 pixel noise. Thus, if the noise power at harmonic £ is {ST^Y from 4c, the upper error 
bar in 4a and 4b will extend from to {T^ + STf^ where Tf is the observed V-W anisotropy 
of each real data point (given by the intersection of the error bars with the zig-zag line) in 
4a and 4b. The rising trend (~ /^) of all three curves towards higher / simply refiects the 
relatively larger pixel noise for smaller angular areas. For / > 200 the real data of 4a and 4b 
rapidly become noise dominated. 
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Fig. 5. — Upper map shows 21 cm data of HVCs with HI column density (A^hi) larger than 
7 X 10^*^ cm~^ (i.e. the greyscale shows A^hi with the outer contour at 7 x 10^^ cm^^). 
Complex C is the cloud in the region / = 90° - 130°, b = 40° - 60°. Complex A is around 
/ = 150°, 6 = 30° - 45°. The Magellanic Stream (MS) and Bridge is at / = 280° - 310°, b < 
-30°. The Leading Arm of the MS, plus some other bright HVCs are at / = 240° - 300°, 
b = 10° - 30°. Lower map gives our estimated V-W temperature excess due to HVCs. Note 
that because the dynamic range of conversion from A'hi to this excess (via free-free emission 
measure EM of A^hii) is not linear (e.g. Putman et al 2003, Hill et al 2009). Our approach 
is to assign 0.5 and 1.0 unit of EM, or 0.15 and 0.3 /xK of V-W excess, to every direction 
with A^Hi > 2 X 10^^ cm~^ and 5 x 10^^ cm~^ respectively. 



